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Abstract 

We importantly amend a certain parenthetical remark made in Part I (arXiv:0806.3294) , to the 
effect that although two-qubit diagonal-entry-paiameteuzed separability functions had been shown 



(arXiv:0704.3723| to clearly conform to a pattern dictated by the "Dyson indices" (/? = 1 [real], 2 
[complex] , 4 [quaternionic] ) of random matrix theory, this did not appear to be the case with regard 
to ei(7en?;a/ue-parameterized separability functions (ESFs). But upon further examination of the 
extensive numerical analyses reported in Part I, we find quite convincing evidence that adherence to 
the Dyson- index pattern does also hold for ESFs, at least as regards the upper half-range ^ < C < 1 
of the maximal concurrence over spectral orbits, C = C(Ai . . . A4) = max{0, Ai — A3 — 2\/A2A4}, 
Ai > A2 > A3 > A4, with the A's being the eigenvalues of associated 4x4 density matrices. To 
be specific, it strongly appears that in this upper half-range, the real two-qubit ESF is simply 
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proportional to (2 — 2C)2 , and its complex counterpart~in conformity to the Dyson-index pattern- 
proportional to the square of the real ESF, that is, (2 — 2C)'^. The previously documented piecewise 
continuous ( "semilinear" ) behavior in the lower half-range < C < I stih appears, however, to 
lack any particular Dyson-index-related interpretation. 
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Part I of this study [T] had been devoted to the question of determining for the generic 
(9-dimensional) real and (15-dimensional) complex two-qubit systems, the nature of cer- 
tain trivariate "eigenvalue-parameterized separability functions" (ESFs). These (metric- 
independent) ESFs, it was argued, could substantially assist in the determination of separa- 
bility probabilities in terms of certain metrics (the Hilbert-Schmidt and Bures being the most 
conspicuous examples). We further investigated in [1] the possibility that these prima facie 
trivariate functions of the eigenvalues Aj (z = 1, . . . 4) of 4x4 density matrices (A4 = 1— Aj), 
were expressible as univariate functions 

sf{X,...X,)=a^^\C{X^...X,)), (1) 

of the maximal concurrence C over spectral orbits [21 sec. VII] [3],|1], 

C(Ai...A4) =max{0,Ai-A3-2v/A^}, Ai > A2 > A3 > A4. (2) 

(At this starting point in our presentation, let us regard /5 in ([T| only as a notational [dummy 
variable] , not calculational device-motivated by Dyson-index conventions-taking the values 
1 [real], 2 [complex], 4 [quaternionic].) 

Our main conclusions in [1] were that-if the reducibility-to-univariance property ([T]) held, 
as our extensive numerical evidence appeared to suggest could be the case-the associated 
real and complex univariate functions both had jumps of approximately 50% magnitude at 
C = ^, as well as a number of additional discontinuities (remarkably coincident in both the 
real and complex cases) in the /ower half-range C G [0, |]. Also, both univariate functions 
appeared to be simply linear between certain of these discontinuities. The upper half-range 
C G [^,l]-in which the univariate functions of C took lesser values-did not command our 
attention in [1], seeming to be of relatively less interest. Our only pertinent observation 
there was that there did not appear to be any discontinuities in this segment. 

Now, in fact, turning our attention more closely to this upper half-range, we readily find 
strong evidence for a very interesting Dyson-index-type phenomenon. If we normalize our 
extensive numerical estimates from p]| of cr^(C) and cr^(C) to both equal 1 at C = |, then 
a joint plot (Fig. [T]) of the latter normalized (complex) function versus the square of the 
former normalized (real) function for C G [|, 1] remarkably shows no perceptible difference 
between the two resulting curves. (The sample [quasi-Monte Carlo] estimate of c^comp/ex^i-j 
is 0.0651586 and that of cr'''^"'(|)) is 0.1803748.) In Fig. |2| we show-on a much finer scale 
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FIG. 1: Joint plot of numerical estimates of j and ^..^^iJ^i'^ for C G [5' 
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than used in Fig. [T]-the actual (very small) differences 
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between them. Of further considerable importance, Fig. |3]is a repetition of Fig. [l| but along 
with the insertion now of the function 



{2 - 2CY = 8(1 - Cf 



(4) 



which we see fits our two estimates very well. Assuming that ^ is the correct form (up to 
the still not exactly-known normalization factor) of ct^(C) over C G 1], we can estimate 
the associated contribution from density matrices corresponding to this half-range to the 
Hilbert-Schmidt and Bures separability probabilities of generic complex two-qubit systems 
to be 0.041568 and 0.0267378, respectively. (The real counterparts of these separability 
probabilities are, then, 0.134611 and 0.104113, respectively.) 
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FIG. 3: The two functions in Fig. [T] along with the additional [very closely-fitting) function 
(2 - 2Cf 



Let us further note that our sample estimate of the ratio 
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0.0651586 
0.18037482 



2.00272 



(5) 



is very close (and possibly theoretically exactly equal) to 2. 

Over < C < i, the range of primary interest in |lj, the estimates of the real and 
complex two-qubit separability functions intersect (near C = 0.1812), and appear to have 
linear segments over the same subintervals [H Figs. 1, 5, 7]. These features, of course, make 
any obvious application of the Dyson-index pattern problematical in this half-range. So, 
the behaviors of the univariate functions a^{C), {/3 = 1 [real], 2 [complex]), over the two 
indicated regimes of C seem to be highly distinct. The point C = ^ clearly serves as a point 
of major behavioral transition, with the lower half-range now appearing perhaps to be the 
more theoretically challenging of the two. An outstanding question would seem to be what 
are the specific values of a'^^°'\^) and o-compiex^ i-j^ which we used as normalization factors in 
our analyses above. The nearness to 2 of the ratio ([s]) may be a helpful guide in this regard. 

Our analyses of two-qubit (img'ona/- entry-parameterized separability functions O EJ [7] 
and eigenvalue-paiameterized separability functions [H |H] have shared a common goal: the 
determination of two-qubit separability volumes and probabilities (in terms of various met- 
rics). As pieces of these formidable objectives begin to be assembled, we can pose a further 
challenge-to find transformations between the two different sets of coordinates used-that 
is, (1) the diagonal entries and (2) the eigenvalues of 4 x 4 density matrices-that will map 
one set of separability functions into the other. The Schur-Horn Theorem, which asserts 
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that the decreasingly-ordered vector of eigenvalues of an Hermitian matrix majorizes the 
decreasingly-ordered vector of its diagonal entries P, chap. 4] (cf. [lOj), would appear to 
be of possible relevance in this regard, particularly since the maximal concurrence C over 
spectral orbits ^ is expressed in terms of the ordered eigenvalues. 
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